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ABSTRACT 

In order to understand what controls the star formation process in luminous 
starburst mergers (e.g., NGC 6240, Arp 220, and so on), we investigate 
observational properties of two samples of high-luminosity starburst galaxies 
mapped in CO(J=1-0) independently using both the Owens Valley Radio 
Observatory (Scoville et al. 1991) and the IRAM interferometer (Downes & 
Solomon 1998). We find that the surface density of far-infrared luminosity, 
S(FIR), is proportional linearly to the H 2 surface mass density, S(H 2 ), for the 
two samples; E(FIR) oc S(H2) 101±0 ' 06 with a correlation coefficient of 0.96. It 
is often considered that S(FIR) provides a good measure of the star formation 
rate per unit area, S(SFR). It is also known that molecular gas is dominated in 
circumnuclear regions in the luminous starburst mergers; i.e., S(gas) ~ S(H 2 ). 
Therefore, the above relationship suggests a star formation law; S(SFR) 
oc S(gas). We suggest that this star formation law favors the gravitational 
instability scenario rather than the cloud-cloud collision one. 

Subject headings: galaxies: starburst - ISM: molecules - stars: formation 
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1. INTRODUCTION 

Since the discovery of ultraluminous infrared galaxies (ULIGs: Soifer et al. 1984; 
Wright, Joseph, & Meikle 1984; see for a review Sanders & Mirabel 1996), much attention 
has been paid to this class of galaxies (hereafter starburst mergers^). Most of these galaxies 
show morphological evidence for galaxy mergers or strong galaxy interaction. There is still 
a controversy on the origin of their huge infrared luminosities; i.e., a pure starburst or a 
hidden central engine of active galactic nuclei (AGN). Recently, growing evidence has been 
accumulated for the relative importance of starbursts in their central regions (e.g., Shaya et 
al. 1994; Genzel et al. 1998; Taniguchi, Trentham, & Shioya 1998 and references therein). 
However, an important question still remains; what controls the star formation process in 
the starburst mergers ? Therefore, it is very important to investigate this problem. 

A number of previous CO observations have shown that the starburst mergers have 
a lot of molecular gas up to ~ 1O 1O M (Young et al. 1986a, 1986b; Sanders et al. 1988; 
Solomon & Sage 1988; Sanders, Scoville, & Soifer 1991). Moreover, the molecular gas is 
usually distributed in the inner 1 kpc region in these mergers [Scoville et al. 1991 (hereafter 
S91); Scoville, Yun, & Bryant 1997; Downes & Solomon 1998 (hereafter DS98)]. Therefore, 
it has been considered that the accumulation of a lot of molecular gas toward the central 
region can lead to the intense starbursts in the mergers. Indeed, numerical simulations 
have shown that galaxy mergers can fuel gas clouds toward the central region of the 
merger products (Barnes & Hernquist 1991; Mihos & Hernquist 1994). Here a question 
arises as how the intense starbursts were triggered in the fueled gas. There may be three 
alternative basic physical processes to cause the intense starbursts; 1) the gravitational 
instability in nuclear gas disks (e.g., Kennicutt 1998 and references therein), 2) efficient 

Tn this Letter, we do not use a more popular terminology, ultraluminous infrared galaxies, 
because our sample studied here contains mergers with L(FIR) < 1O 12 L . 
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cloud-cloud collision driven by violent gas motion in their central regions (e.g., S91), and 
3) triggering through the dynamical disturbance driven by supermassive binary black holes 
(Taniguchi & Wada 1996). In order to investigate which process is most important in the 
starburst mergers, we analyze two data sets of high-luminosity starburst mergers mapped 
in CO(J=1-0) using the Owens Valley Radio Observatory (OVRO) interferometer (S91) 
and the IRAM one (DS98). 



2. DATA 

The first data set is taken from S91 who presented a summary of the OVRO CO( J =1-0) 
mapping of 18 high-luminosity infrared galaxies. Almost all the galaxies show evidence 
either for galaxy mergers or for strong galaxy interaction. This sample covers a wide range 
of far-infrared (FIR) luminosities from ~ 1O 1O L to ~ 1O 12 L . The second data set is 
taken from DS98 who performed detailed CO( J=l-0) and CO( J=2-l) observations of ten 
ULIGs. Though four ULIGs among this sample are also included in the OVRO sample (Arp 
220, Mrk 231, VII Zw 31, and IRAS 17208-0014), we deal with this IRAM sample as an 
independent sample appreciating the two independent CO mapping programs using the two 
different radio interferometers. Although these two samples are not statistically complete 
ones, they are highly useful in investigating a star formation law in starburst mergers. In 
Table 1, we give a summary of observational properties of the individual galaxies of the 
two samples. We use the distances of galaxies estimated using both recession velocities 
corrected to the Galactic Standard of Rest (GSR) following the manner described in de 
Vaucouleurs et al. (1991) and a Hubble constant H = 75 km s -1 Mpc -1 . The heliocentric 
velocities are taken from de Vaucouleurs et al. (1991), S91, and DS98 (see Table 1). 

Molecular gas masses for all the galaxies in the two samples are derived using a 
conversion factor, a = Ah 2 //(CO) = 3 x 10 20 cm -2 (K km s -1 ) -1 where Nu 2 is the H 2 
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column density and /(CO) is the integrated intensity of the CO(J=1-0) emission (S91). 
The FIR data are compiled from the IRAS Faint Source Catalog (Moshir et al. 1992). The 
FIR luminosities are estimated using L(FIR)= 4vr J D 2 1.26 x Kr n [2.58 x 5(60) + 5(100)] 
(ergs s _1 ) where 5(60) and 5(100) are the IRAS 60 /im and 100 /xm fluxes in units of Jy 
and D is the distance of galaxies in units of cm (Helou, Soifer, & Rowan- Robinson 1985). 
Finally L(FIR) is given in units of L Q . In later analysis, in addition to L(FIR) and Mh 2 , 
we also use their surface densities; S(FIR) (L Q pc~ 2 ) and £(H 2 ) (M Q pc~ 2 ), respectively. 
These densities for each galaxy are normalized by a surface A = nr 2 where r is the radius 
in units of pc given in S91 and DS98. Note that for the four ULIGs commonly observed by 
both S91 and DS98, the CO sizes in DS98 are systematically smaller than those in S91 (see 
Table 1). We use the CO sizes given in DS98 for these four ULIGs in the later analysis. 



3. RESULTS 

In Figure 1, we show diagrams of log L(FIR) (L ), log S(FIR) (L pc~ 2 ), and log 
L(FIR)/M H2 (L g /M q ) as functions of log S(H 2 ) (M pc~ 2 ) and log 5(60)/5(100) for the 
two samples of S91 and DS98. Among the six relationships shown in Figure 1, we find the 
following two significant correlations for both the samples; a) log S(FIR) vs. log X(H 2 ), 
and b) log L(FIR)/M dust vs. log 5(60)/5(100). 

We obtain the nearly linear relationship between S(FIR) and S(H 2 ); S(FIR) 
oc S(H 2 ) L02±a06 for the S91 sample, E(FIR) oc £(H 2 ) L16±(U3 for the DS98 sample, and 
S(FIR) oc S(H 2 ) L01±0 - 06 for the total (S91+DS98) sample. Note that the new data in DS98 
are adopted for the four common ULIGs in the total sample. Therefore, the total number 
of ULIGs in the total sample is 24. The correlation coefficient for the total sample is quite 
high, R ~ 0.96. Since S(FIR) can be regarded as the star formation rate (SFR) per unit 
surface area (e.g., Kennicutt 1998), i.e., S(SFR), this correlation implies that the SFR is 
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controlled locally by the surface gas density. On the other hand, since there is no significant 
correlation between L(FIR) and S(H 2 ), the total FIR luminosity is governed not by the 
local gas density but by the total gaseous content accumulated in the central region of 
starburst mergers. 

The correlation between L(FIR)/M H2 and 5(60) /S( 100) has been already found in a 
number of studies (e.g., Young et al. 1986a). The correlation found here is not as tight 
as that for S(FIR) vs. log E(H 2 ). The data points of both NGC 4038/4039 and IRAS 
00057+4021 make the correlation worse though we do not understand why these data are 
far from the correlation for the remaining galaxies. Since the L(FIR)/M H2 ratio provides 
a measure of star formation efficiency (Young et al. 1986a, 1986b; S91), this correlation 
suggests that the higher star formation efficiency leads to the production of more heating 
photons per dust grain, resulting in the higher dust temperature. 

Although S91 obtained a significant correlation between log L(FIR)/M H2 vs. log S(H 2 ) 
for 14 high-luminosity galaxies, we do not find this correlation for all the galaxies in S91 
and for the galaxies in DS98. The little correlation between L(FIR)/M H2 and S(H 2 ) means 
that the star formation efficiency is not controlled by the surface gas density. 

In Figure 2, we show diagrams of S(FIR) and S(H 2 ) as a function of A in logarithmic 
scales. It is shown that the two kinds of surface densities are almost inversely proportional 
to the surface, A; 1) S(FIR) cx A a : a = -0.77±0.15 for the S91 sample, -1.05±0.15 for the 
DS98 sample, and -0.85 ± 0.13 for the total sample, and 2) £(H 2 ) oc A b : b = -0.74 ± 0.13 
for the S91 sample, -0.90 ± 0.09 for the DS98 sample, and -0.82 ± 0.12 for the total 
sample. The correlation coefficients for the total sample are —0.80 and —0.82, respectively. 
Therefore, one important thing in the luminous starburst mergers is how the gas clouds are 
accumulated in smaller regions. We suggest that more concentrated gaseous systems tend 
to experience more intense starbursts such that S(SFR) oc E(gas). 
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4. DISCUSSION 

We consider the physical meaning of the linear relationship between E(FIR) and X(H 2 ). 
As mentioned before, this relationship suggests the star formation law of S(SFR) oc £(H 2 ). 
This kind of relation has been discussed by many authors in terms of so-called Schmidt law 
(Schmidt 1959); i.e., the star formation law remains the simple gas-density power law; i.e., 
S(SFR) oc ^(gas)^. In the derivation of the Schmidt law, it is recommended to use the 
total (i.e., H I + H 2 ) gas mass density. However, since the nuclear gas in the high-luminosity 
starburst mergers is often dominated by molecular gas (Mirabel & Sanders 1989; see also 
Sanders & Mirabel 1996), we can neglect the contribution of H I gas (Kennicutt 1998); i.e., 
E(gas) ~ S(H 2 ). Therefore, the linear correlation between S(FIR) and S(H 2 ) suggests that 
the star formation in the luminous starburst mergers can be described by a Schmidt law 
with N ~ 1. 

For disks of normal galaxies, Kennicutt (1989) derived a relation of N ~ 1.3 ± 0.3 (see 
also, Buat, Deharveng, & Donas 1989; Buat 1992; Bosseli 1994; Deharveng et al. 1994). 
Recently, Kennicutt (1998) derived a star formation law for a sample of 36 infrared-selected 
starburst galaxies; N ~ 1.40 ±0.13. This Schmidt law is much steeper than that derived in 
this study. This seems as due to the difference in sampling between Kennicutt (1998) and 
ours. In his sample, a number of less luminous starburst galaxies (e.g., IC 342, NGC 6946, 
and so on) are included. Although these galaxies experience the moderate star forming 
activity, they are basically isolated ordinary spiral galaxies. On the other hand, our sample 
comprises mostly major mergers. In fact, the steeper Schmidt law derived by Kennicutt 
(1998) appears to be biased by the inclusion of the less- luminous, non-merger galaxies. 

As mentioned in Section 1, there are three alternative mechanisms to initiate starbursts 
in merging galaxies; 1) the gravitational instability in nuclear gas disks (e.g., Kennicutt 
1998 and references therein), 2) efficient cloud-cloud collision driven by violent gas motion 
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in their central regions (e.g., Young et al. 1986a; S91), and 3) triggering through the 
dynamical disturbance driven by a supermassive binary (Taniguchi & Wada 1996). A 
Schmidt law with N ~ 2 can be achieved by the star formation driven by cloud-cloud 
collisions (Scoville, Sanders, & Clemens 1986). On the other hand, a Schmidt law with 
AT ~ 1 - 1.5 implies the star formation triggered by the gravitational instability of gaseous 
disks (Kennicutt 1989, 1998; Larson 1988; Elmegreen 1994). The Schmidt law derived 
here (i.e., N ~ 1) favors the gravitational instability scenario rather than the cloud-cloud 
collision model. Finally, we cannot rule out the third possibility because this model also 
adopts the gravitational instability formalism. 

Recently, Taniguchi, Trentham, & Shioya (1998) proposed a starburst-driven starburst 
mechanism in order to understand the formation of blue super star clusters in ULIGs such 
as Arp 220 (see also Shioya, Taniguchi, & Trentham 1998). In such a mechanism the 
infalling dense gas disk is unstable gravitationally and collapses to form massive gaseous 
clumps. These clumps are exposed to the external high pressure from the superwind0 driven 
by the ongoing starburst at the very center of the galaxy. This external pressure leads to 
further collapse of the clumps which in turn leads to massive star formation in them. Our 
result is consistent with this scenario. 

Starbursts in major mergers have been also studied by numerical simulations (e.g., 
Mihos & Hernquist 1994 and references therein). In most the simulations with an SPH 
code, the star formation law is often assumed to follow a Schmidt law. For example, Mihos 
& Hernquist (1994) adopted the Schmidt law with iV = 1.5. However, our result shows that 
the Schmidt law with iV ~ 1 describes the star formation in the actual high-luminosity 
major mergers. Therefore, we recommend that future numerical simulations for gas-rich 

2 A blast wave driven by a collective effect of a large number of supernovae in the very 
core of the galaxy (e.g., Heckman, Armus, & Miley 1990; Heckman et al. 1996). 
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major mergers will use the Schmidt law with N ~ 1. 
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Table I. Properties of the sample of galaxies 
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Fig. 1.— Diagrams of log L(FLR), log S(FIR), and log L(FIR)/M H2 as functions of log 
S(H 2 ) and log S(60)/S(100). The OVRO sample is shown by filled circles while the IRAM 
one is shown by open circles. The filled squares for the four ULIGs observed by S91 (Arp 
220, Mrk 231, VII Zw 31, and IRAS 17208-0014) show the data points when we adopt the 
sizes obtained by DS98. However, note that the data points of the four ULIGs in the log 
S(FIR) vs. log S(60)/S(100) are the same as those for DS98. 

Fig. 2. — Diagrams of log S(FIR) and log S(H 2 ) as a function of the surface A. The symbols 
have the same meanings as those in Figure 1. 
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Table 1 — Continued 
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Note. - (1) Name, (2) Distance (Mpc); heliocentric velocities are taken from (b) 

de Vaucouleurs et al. (1991), (c) S91, and (d) DS98, (3) Radius (kpc), (4) S(60) (Jy), 
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L(FIR)/M(H 2 ) (L /M ), (10) log S(FIR) (L /pc 2 ), and (11) log S(H 2 ) (M /pc 2 ) 
a I17208 = IRAS 17208-0014; 110173 = IRAS 10173+0828; Zw049 = Zw049.057; 

100057 = IRAS 00057+4021; 102483 = IRAS 02483+4302; 110565 = IRAS 

10565+2448; 123365 = IRAS 23365+3604. 
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